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We show how a ground state trial wavefunction of a Fermi liquid can be systematically improved 
introducing a sequence of renormalized coordinates through an iterative backflow transformation. 

We apply this scheme to calculate the ground state energy of liquid ^He in two dimensions at freezing 
density using variational and fixed-node diffusion Monte Carlo. Comparing with exact transient 
estimate results for systems with small number of particles, we find that variance extrapolations 
provide accurate results for the true ground state together with stringent lower bounds. For larger 
systems these bounds can in turn be used to quantify the systematic bias of fixed-node calculations. 

These wave functions are size consistent and the scaling of their computational complexity with the 
number of particles is the same as for standard backflow wave functions. 

PACS numbers: PACS: 


I. INTRODUCTION 

To overcome the fermion-sign problem, many fermion 
Quantum Monte Carlo (QMC) calculations rely on the 
fixed-node (FN) approximation where the nodes of a trial 
wavefunction, ipT, are imposed as a boundary condition 
on the many-body Schrddinger equation which can then 
be solved by projector Monte Carlo methodi^^. Since the 
nodal surfaces of the exact ground state wavefunction are 
in general unknown, the energies of FN calculations do 
not converge to the exact ground state energy but remain 
above them by an unknown amount. Although methods 
which do not rely on the FN approximation have been 
developecPlf^, they are in general limited to small systems 
as their computational cost grows exponentially with sys¬ 
tem size. Therefore, FN-QMC calculations still provide 
the most accurate values of ground state properties of 
extended fermion systems. 

Modification of the nodes of a many-fermion wave 
function to explicitly include correlations remains a 
formidable task. Slater determinants based on Backflow 
(BF) coordinates present one possibilitji^^, and back- 
flow wave functions have been routinely used over the 
last years in QMC calculations of the electron-gaJ^i®^ 
and liquid Generalization of the backflow wave 

function to include three body correlations was shown to 
be necessary to stabilize the unpolarized phase of liquid 
^He against spin-polarizatiorfi^. 

Here we propose new correlated trial wave functions 
based on iterative backflow transformations and use them 
to study liquid ^He in two dimensions. We show that this 
new class of trial wave functions systematically lower the 
energy and its variance. Our results illustrate the pos¬ 
sibility to extrapolate variational (VMC) and FN diffu¬ 
sion (DMC) Monte Carlo calculations to zero variance 


to approach very closely the exact ground state energy. 
Since their evaluation remains of similar complexity and 
scaling with increasing system size as the usual backflow 
wave function, their use is not limited to small systems. 
We explicitly demonstrate the size consistency of our new 
trial wave functions and discuss the possibility to obtain 
lower bounds to the ground state energy. 


II. ITERATED RENORMALIZATION OF WAVE 
FUNCTION 

Let us start by considering the standard Slater-Jastrow 
type trial wave function with backflow, 

= det(/ifc(q4R])e“^[^l (1) 

Antisymmetry is ensured by the Slater determinant of 
single particle orbitals, (/)fc(r), k = where, in¬ 
stead of the bare coordinates r^, i = many- 

body backflow coordinates, q^, are used as arguments. 
Both backflow coordinates, Q = (qi,..., qAr), and the 
symmetric Jastrow potential, U, depend explicitly on all 
coordinates, R = (ri, r 2 ,..., rjv), as indicated. In the 
standard form, U = + SiGi(R) • Gj(R) 

with G* = Ej (u-rj)C(r„ ), q* = U-k Ej (u - rj)7?(ry), 
and rij = jr^ — r^j. The radial functions u, ^ and ij 
can be parametrized and optimized by minimization of 
the variational energy. Generalizations to include higher 
correlations into both backflow and Jastrow potentials 
are possibl^i^ but will not be considered here. 

Once the backflow and Jastrow potentials have been 
determined, different occupations of the orbitals inside 
the Slater determinant of Eq. 0 can be used to ap¬ 
proximate also low-lying excited states of the systems^. 
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in close analogy to Landau’s Fermi liquid de script ion. 
As in the correlated basis functions approac H^^ I ^^ I, let 
us consider the effective Hamiltonian within these non- 
orthogonal basis states. For a Fermi liquid, we expect 
non-diagonal matrix elements of the effective Hamilto¬ 
nian to be strongly suppressed compared to those of the 
bare plane-wave states. However, instead of diagonaliz¬ 
ing the effective Hamiltonian, let us search again for a 
trial wave function to represent the ground state of the 
effective Hamiltonian. Assuming a smoothly varying ef¬ 
fective interaction, we may again consider to represent it 
as a backflow wave function, However, this time, 

the new backflow coordinates, and the new Jastrow 
potential, are built upon the old backflow coordi¬ 
nates, with = Q. Thus 

we are naturally led to an iterative renormalization pro¬ 
cedure 

^ = det (/>fe(qf(2) 

with a renormalized Jastrow potential 


III. TWO DIMENSIONAL LIQUID ’^HE AT 
FREEZING DENSITY 

In order to illustrate the accuracy of the renormaliza¬ 
tion procedure to describe the ground state wave function 
of highly correlated Fermi liquids, we perform calcula¬ 
tions for the ground state energy of liquid ^He in two 
dimensions at a density p = 0.060A“^, near freezinj^. 
We compare VMC and fixed node DMC energies to ex¬ 
act results obtained by the nominally exact transient es¬ 
timate (TE) method of Ref. [3l for systems of IV = 26 
{N = 29) unpolarized (polariz^ ^He atoms interact¬ 
ing with the HFDHE2 potentiaP^l. Furthermore we test 
the size-consistency of our trial functions, comparing the 
gain in variational energy obtained by the renormaliza¬ 
tion procedure for the unpolarized system at two different 
sizes, N = 26 and N = 58. The results are collected in 
Table H 

Every iteration introduces three new potentials (for 
backflow, two- and three-body Jastrow function), each of 
which, generically indicated here as /(r), is parametrized 
in the form 


= (,!<>->) 

/3<a i<j 


( 3 ) 


and renormalized backflow coordinates 

q'“’ = r. + E y.'*” (4) 

l3<ct 

= E (q‘“-' - qf-’) ■)'"> (?‘r‘>) 

3¥=i 

(in Eqs. ^ and Q, stands for R). At each it¬ 

eration new potentials parametrizing the additional Jas¬ 
trow and backflow functions are introduced, and all the 
potentials and with /? < a, have to be 

optimized. 

In the appendix we show how the evaluation of the 
renormalized wave functions and their derivatives needed 
to calculate the local energy can be efficiently imple¬ 
mented with a number of operations proportional to N^. 
Thus, the overall cost of calculation is not dramatically 
altered compared to the usual (zeroth order) backflow 
wave function. For a system of = 26 particles we find 
that the CPU time to move all the particles and calculate 
the local energy with iterated backflow of order 1 to 4 is 
a factor 5, 9, 13, 17 larger than that of the zeroth order, 
respectively; furthermore, for N = 58, fourth order back- 
flow takes 12.5 times longer than for N = 26, close to the 
scaling. The corresponding figures for the efficiency 
of the calculation of the energy are even more favorable, 
because the variance is lower for improved wave func¬ 
tions. 


fir) 


CO 

1 

0 

En=l ^ + 06/?’“" 

[o 



if r < rc 
if r>rc- 


( 5 ) 

For the backflow and three-body Jastrow potentials we 
set rc = TA and drop the McMillan term {qq = 0), while 
for the two-body Jastrow potential we choose a cutoff 
value rc close to half the side of the simulation box. In 
Fig. 0 we show the optimized potentials and 

of the 4*^^ wave function for a system with N = 26 
and C = 0. The backflow coordinate transformations 
across different iterations implicitly build up many-body 
correlations at all orders, so that eventually not all of 
the optimized potentials have an obvious physical inter¬ 
pretation: for instance the pair distribution functions g 
of the bare coordinates and of the renormalized coordi¬ 
nates at subsequent iteration levels, shown in Fig. fea¬ 
ture increasingly wide correlation holes and high peaks 
for increasing level, despite the two-body potentials 
turning from repulsive for a = 0 to attractive for a = 4. 
Note that all the g's feature the structure of simple liq¬ 
uids (albeit with increasingly classical character), which 
supports the heuristic derivation given in Section [IT} each 
iteration essentially renormalizes the Slater Jastrow wave 
function without qualitative changes. 

With the choice of Eq. ([^, the renormalization pro¬ 
cedure requires 17 variational parameters per level, and 
the corresponding optimization procedure (carried out by 
correlated sampling^^ in this work) becomes rather de¬ 
manding. Therefore we have tried two simpler iterative 
schemes, one in which no renormalized Jastrow is present, 
and one in which only the new potentials added at the a- 
th iteration are optimized, leaving the other unchanged 
from previous iterations. However, these simpler options 
lead to higher values in energy, both for VMC and DMC. 
We have also considered an improved wave function with 
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N=26, C = 0 

N=58, C = 0 

N=29, C = 1 


Et/N 

a^/N 

A 

Edmc/N 

Et/N 

M/N 

A 

Edmc/N 

Et/N 

M/N 

A 

Edmc/N 

PW 

3.011(1) 

28.29 


2.419(2) 

2.900(1) 

28.07 


2.373(2) 

2.5831(6) 

7.51 


2.402(1) 

BF(°1 

2.688(1) 

13.05 

0.323 

2.353(2) 

2.584(1) 

13.34 

0.316 

2.283(2) 

2.5133(5) 

5.34 

0.070 

2.4005(6) 

BF^il 

2.471(1) 

4.58 

0.540 

2.336(2) 

2.356(2) 

4.93 

0.544 


2.4383(3) 

2.20 

0.145 

2.3918(5) 

BF(2) 

2.4258(8) 

2.86 

0.585 

2.3284(9) 

2.313(2) 

3.25 

0.587 


2.4193(3) 

1.54 

0.164 

2.3877(4) 

BF'^1 

2.4049(9) 

2.47 

0.606 

2.3223(4) 

2.297(2) 

2.67 

0.603 


2.4136(2) 

1.36 

0.170 

2.387(1) 

bf('‘i 

2.400(1) 

2.29 

0.611 

2.323(1) 

2.292(1) 

2.49 

0.608 

2.232(1) 

2.4109(7) 

1.25 

0.173 

2.3869(5) 

VMCe.t 


2.338(5) 



2.217(2) 



2.384(6) 


DMCerit 


2.317(3) 



2.216(3) 



2.379(1) 


LBerct 


2.275(14) 



2.149(12) 



2.390(26) 


TE 

2.307(7) 


2.375(3) 


TABLE I. Ground-state energy per particle, in K, of liquid ^He in two dimensions at p = O.O 6 OA obtained with Variational 
(Et/N) and fixed node Diffusion Monte Carlo {Edmc/N) using different types of trial wave functions: Slater-Jastrow wave 
function without backflow (PW), and with Q-times iterated backflow is the spin polarization and N is the number 

of particles. A is the gain in VMC energy per particle relative to the PW value, and is the variance of the VMC total 
energy. TE indicates unbiased results calculated with the transient estimate method of Ref. [S] VMCext, DMCea;t and EB^xt 
are the extrapolations to zero variance of Et/N, Edmc/N and of the lower bound {Et — '/(E)/N, respectively. Statistical 
uncertainties on the last digit(s) are given in parentheses. All values are given for periodic boundary conditions (P point) 
without tail corrections^. 


different backflow potentials for parallel and antiparallel 
spins. The gain in energy is ^ 10 mK in VMC, but hardly 
visible in DMC (< 1 mK) beyond the second backflow 
iteration. Finally, we have tested the accuracy of using 
potentials optimized for TV = 26 to perform simulations 
with N = 58 particles: at the fourth backflow iteration, 
the DMC energy is higher by a non negligible amount, 
5±2 mK. All these results are listed in Table 0 

IV. ZERO-VARIANCE EXTRAPOLATION AND 
LOWER BOUNDS 

Our VMC and FN-DMC results for the energy expec¬ 
tation values Ex = of the different (nor¬ 

malized) wave functions, provide strict upper bounds 
for the true ground state energy, Eq < Ex, where the 
subscript X stands for T or DMC as appropriate and 
^DMC is the FN ground state. Within VMC, we have 
also access to the variance of the energy in the trial 
state, — Et)'^\'^t)- As the variance ap¬ 

proaches zero for any exact eigenstate, its value for a 
given trial wave function can be used to quantify the 
distance to the closest eigenfunction. Under the assump¬ 
tion that the trial energy is closer to the ground state 
energy than to any of the other eigenstates, the inequal¬ 
ity tr^ > {Eq — Et)^ leads to a lower bound for the ground 
state energji^ 


Eo>Et- VctA (6) 

In the following, we will use the information on the 
variance obtained by VMC to extrapolate to the exact 
ground state energy. 

Let us first analyse in more detail how the trial wave 
function approaches the ground state wave function. Ex¬ 
panding our trial wave function in the exact eigenstates. 


\Ej), of energy Ej, we have where Cj 

are the expansion coefficients, with = 1 assum¬ 

ing normalized states. We can now write 


Et = Eq XtCt (7) 

with Ct = P and At = J2iiEi-Eo)cf/C t > A, 

where A = Ei — Eq denotes the energy gap between 
ground and first excited state of the system. Similarly, 
we obtain for the variance 

(T^ = A'^tCt — {AtCt) ( 8 ) 


where A'^t = “ ^o)^cf/C'r > Af,. 

Using Eqs. Q and Q we have 


Et — a / = Eq — AtCt 


I A^t 
AICt 


- 1-1 


( 9 ) 


and we see that the expression for the lower bound, 
Eq. 1^, remains valid for Ct < A'^t/‘^A\, or Et — Eq < 
A'^t/'^At- Note that this condition is less stringent than 
the assumption that Et is closer to the ground state en¬ 
ergy than to any of the other eigenstates used previously. 

To go further, let us assume that the trial wave func¬ 
tion has a significant overlap only with the ground state 
wave function, whereas the components of excited staes, 
a with f > 0, are broadly distributed. We expect this 
assumption to be reasonably satisfied for extended sys¬ 
tems, where the excited states approach a continuum in 
the thermodynamic limit. Improving the wave function 
via our iterative renormalization, the excited state con¬ 
tributions decrease almost uniformly, such that Ct —t 0 
whereas At and A^t remain roughly constant. In this 
case we can neglect terms of order C^ in Eq. ([^ and 
insert it in Eq. Q, and obtain 

Et = Eq + Aa'^, for cr^ —)• 0 


( 10 ) 
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iV = 26, C = 0 





A = 58, C = 0 




Et/N 



Edmc/N 


Et/N 

Edmc/N 


0 

I 

II 

III 

0 

I 

II 

III 

0 

IV 

0 

IV 

PW 

BpC) 

Bpl^l 

Bpl^) 

2.471(1) 

2.4258(8) 

2.4049(9) 

2.599(2) 

2.585(2) 

2.584(2) 

2.515(1) 

2.480(1) 

2.472(1) 

2.461(1) 

2.413(1) 

2.398(1) 

2.336(2) 

2.3284(9) 

2.3223(4) 

2.337(2) 

2.335(2) 

2.335(2) 

2.337(2) 

2.332(1) 

2.326(1) 

2.3256(9) 

2.3215(4) 

2.900(1) 

2.584(1) 

2.909(2) 

2.592(2) 

2.289(2) 

2.288(1) 

BpU) 

2.400(1) 

2.580(2) 

2.470(1) 

2.390(2) 

2.323(1) 

2.331(1) 

2.325(3) 

2.324(1) 

2.292(1) 

2.298(2) 

2.232(1) 

2.237(1) 


TABLE II. Some of the energies of Table |I] compared to the corresponding values obtained with downgraded or upgraded 
wave functions. Entries 0: energies from Table [T] entries I: downgraded wave functions with omitted Jastrow factors of 
the quasi-coordinates; entries II: downgraded wave functions with Jastrow and backflow potentials from previous iterations 
not repotimized; entries III: upgraded wave functions with different like-spin and unlike-spin backflow potentials; entries IV: 
downgraded wave functions for N = 58 with Jastrow and backflow potentials optimized for N — 26. 


with A = Therefore, with good enough trial 

functions, we expect that a (nearly) linear extrapola¬ 
tion of the variational energy to zero variance closely 
approaches the exact ground state energy, with the coef¬ 
ficient of the linear term providing a numerical estimate 
of the validity of the lower bound of Eq. §, i.e. 

Et-Eo< 1/(2A). (11) 

The lower bound, in turn, can be made stricter by ex¬ 
trapolation to zero variance of Ej- — with a lead¬ 
ing square-root term. These variance extrapolations are 
shown in Figs Q and 0 and listed in Table |!) 

All of the above extrapolations are valid for Ct —t 
0. Using a general estimate for the overlap of the trial 
wave function with the ground stat^^, Gt = 1 — Cq > 
1 — exp[—(F t — £'o)^/2cr^] « 1 — exp[—A^cr^], we can 
a-posteriori check the consistency of the energy versus 
variance extrapolation, see Fig. 

Finally, one would like to use variance extrapolation 
with the fixed-node energies to obtain even better re¬ 
sults. However, within DMC, the variance dmc\{E[ — 
Edmc)‘^\^DM c) is zero inside any nodal pockeP^ and 
cannot be used anymore as a measure of the quality of 
the wave function. The most natural assumption is then 
to postulate that the variance calculated in VMC is 
a good measure of the quality of the wave function in 
DMC as well. This allows us to use the same extrapola¬ 
tion for DMC energies as in the case of VMC, as shown 
in Figs. @, 0 and 0, but without obtaining a lower 
bound. The DMC energy extrapolated to zero variance, 
listed in Table [ij happens to differ from the TE value by 
just the combined error bar, 10 mK for C = 0 and 4 mK 
for C = 1. 

V. CONCLUSIONS 

In this paper we have introduced new, highly corre¬ 
lated wave functions for accurate descriptions of normal 
Fermi liquids based on generalized backflow coordinates 
which are iteratively improved. For liquid ^He at freez¬ 
ing density, the energy gain of these wave functions at 


the 4th iteration compared to the usual backflow trial 
wave function (0 iterations) is about 290mK within VMC 
and 30mK for FN-DMC. More important, we have shown 
that the true ground state energy can be obtained by 
variance extrapolation with intrinsic a-posteriori checks 
of the consistency and validity of the extrapolation. For 
small number of atoms, N ^ 26, we have shown that 
the obtained results are in agreement with unbiased cal¬ 
culations using transient estimates, but variance extrap¬ 
olation can be used also to quantify the fixed-node er¬ 
ror of larger systems. For systems with N = 58 atoms, 
the fixed-node error of our best wave function is around 
20mK. 

Thus, apart from significant VMC and FN-DMC en¬ 
ergy gains, the iterative backflow renormalization pro¬ 
cedure also leads to a general strategy to quantify the 
fixed-node error of the calculations. In combination with 
finite-size extrapolations based on the an alytic al informa¬ 
tions contained in the trial wave functioiisSHU^ the meth¬ 
ods presented in the paper provide an important step 
towards the control of the accuracy of QMC calculations 
suffering from a Fermion sign problem. 
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Appendix A: Computational details 

Let us suppose that gf, are backflow coordinates i = 
1,... N and a = 1,... d, where d is the spatial dimen¬ 
sion, and we have already computed the following partial 
derivatives 

Qtf = (Al) 

= (A2) 

ia 
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We will further need 




(A3) 


where u denotes the function and u' (it") its first (second) 
derivative. The gradient of the Jastrow factor can then 
be calculated by the chain rule 


which is already needed for computation of the local en¬ 
ergy of the Slater determinant using orbitals based on 
the above backflow coordinateJi^^, which we will shortly 
remind. 


1. Backflow determinant 

The gradient and the laplacian of a determinant, D = 
detipki, with backflow coordinates in the orbitals, (p^i = 
can be calculated as follows 


313 


Q.0 

‘i’j 


(A4) 


Alog D = Y,K^Qf + Y. 

ioL ick 






ia.p L ^ 

j'i' 




(A5) 




i/3 


m^l 


n 

Qlm 


and 


AU = J2vfQf 

W 


EEw'i 

l^m /37 


Im 


_ qM 
VZZ V/m 


(A9) 


(AlO) 


with 


- 






u'{qim) 


Qlm J Qlm Qlm 

(All) 

We see that the overall cost of the quasi-particle Jastrow 
factor and its derivatives needed for the local energ y is o f 


order of N^, needed to build the matrix Qf^, Eq. (A31. 


Since this matrix is already needed in the calculation of 
usual backflow wave functiorfS^, the iterated Jastrow does 
not lead to a significant slow down compared to the usual 
backflow. 


where 


„ _ dipki a/3 _ d (fiki pa _ rr , a 


dq° 


and Vik is the inverse of the backflow matrix 
1 dD 

k 


V,k = 


D dipk, 


^ ^ ^ikPkj — 


'3 

(A6) 

(A7) 


The computational complexity is of order for the in¬ 
version of the orbital matrix, V)fc, as well as for the com- 

putation of the matrices and . 

Note that this part of the calculations does not de¬ 
pend on the specific form of the backflow coordinates. 
The computation of the gradient and laplacian of the 
local energy based on Eq (A4) and Eq (A5) do only de¬ 


pend on the actual values of the orbital matrix, pki, and 
its partial derivatives, Eqs. (A6|, and on the gra dient 
and laplacian of the backflow coordinates, Eq (All and 
Eq. (A2). Therefore, Eq. (|A4|) and Eq. (A5| can be still 


used to calculate the local energy of determinants con¬ 
taining iterated backflow coordinates, as long as their 
derivatives are provided in the form of Eqs (Al), (A2|, 


and (A31. 


2. Iterated Jastrow correlations 

We can now build a Jastrow factor based on the dis¬ 
tances between two quasi-particles. 




(AS) 


3. Iterated backflow coordinates 


We now construct new backflow coordinates 

y? = 


(A12) 


where rj is the corresponding potential. In order to cal¬ 
culate the local energy for backflow orbitals in the Slater 
determinant based on ?/“, we need the following deriva¬ 
tives 

y^f = = AyT (A13) 

In order to calculate them, we will use the chain rule, 
based on the following partial derivatives 

EEl - . V iy“'^ - 17“^ 

Q„a — ytn 

T) 


dqfdq^ 

where 


— ‘’ijk2_^ykn ^jktlki “yVjfe ^i-^ykj 


.a/3 _ V'{qij)„a 


V^[ = 


“OcB'y 

y^3 = 


q^J 


yfjQii + ^(9b)^a/3 


v"iQ^3 - 


v'iQij) 


qij 


q“-oE qT 


-b 


q^J 


qfjSp^ + qf ja-y -b qJJcfi 


l<.m 


and we have used that 7 )“^ = 1 /“^"’' = 0. 
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The final derivatives needed, Eqs. (A13), can then be 
written as 


y^f-T.yln[Q7-Q 


a-yl 
in J 


nj3 


-T.ytn" 

noL^ 




(A14) 


Again, these operations can be done in order of com¬ 
putations. 


4. Iterated n-body correlations 

Above we have explicitly shown how to calculate gradi¬ 
ent and laplacian of a scalar two-body Jastrow potential 
and of quasi-particle coordinates constructed from back- 
flow coordinates. The structure of our three-body corre¬ 
lation in Eq. (|^ is actually a scalar product between two 
vectors with identical structure as the backflow coordi¬ 
nates. Gradient and laplacian of the three body term can 


therefore be calculated from those of the vectors using 
the chain rule. Generalizations to build iterated many- 
body Jastrow and backflow coordinates based on quasi¬ 
particle tensorare straightforward and do not increase 
the complexity of the calculation. 


5. Higher order iterations 

At the zeroth order of iteration, the backflow coordi¬ 
nates, qi, are symmetric functions of the bare coordi¬ 
nates. Higher order iterations of the backflow are built 
from symmetric expressions based on the previous back- 
flow coordinates, such that the overall wave function re¬ 
mains antisymmetric. Above, we have explicitly shown 
how to calculate the gradient and the laplacian of the 
quasi-particle coordinates of the hrst backflow iteration, 
= qi -I- Yi without increasing the overall complexity 
of the calculation. These are the only additional infor¬ 
mation needed to calculate the local energy of the first 
iterated backflow determinant, and from the structure it 
is clear that this procedure can be iterated to higher order 
without increasing the complexity of the calculations. 
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FIG. 1. Optimized potentials of the trial function for 
N = 26, = 0. The lines are broken where the pair dis¬ 

tribution functions of the relevant (quasi)coordinates become 
negligibly small, g{r) < 10“^ (see Fig. [^. 
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r (A) 


FIG. 2. Pair correlation functions calcnlated in a VMC simu¬ 
lation with the trial function using the bare coordinates 

{ri} (solid line), and the renormalized backflow coordinates 
}, and (dotted, dash-dotted and dashed 

lines, respectively). The statistical noise reaches its maxi¬ 
mum value ~ 0.003 at the highest peak. 
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FIG. 3. The VMC and DMC energies per particle of A'" = 26 
nnpolarized ^He atoms in two dimensions as a function of the 
variance jN. Each point corresponds to a different trial 
function (PW and BF^“^ with a = 0 to 4, from higher to 
lower variance). The TE energy has been subtracted. Their 
dependence is nearly linear, and their extrapolations to zero 
variance, the entries VMCe3;t and DMCe3;t of Table|I] are very 
close to the exact result. We further show the energy lower 
bound [Et — whose extrapolation to zero variance, 

the entry LBe^t of Table [I] is also very close to the exact re¬ 
sult. The dashed line is a rough estimate of the first excited 
state (the difference between the two slowest exponential de¬ 
cay constants in the fermionic signal of the TE procedur^. 
It shows that the condition for the validity of the energy lower 
bound (see text) are met by the iterated backflow trial func¬ 
tions. The dotted line is the alternate estimate 1/{2AN) of 
Eq. ^ for the validity of the lower bound. 



FIG. 4. Same as Fig. for N = 29 spin-polarized ^He atoms. 
Both the estimate of the first excited state, ~ 0.6 K, and the 
value of 1/{2AN), 0.772 K, are off scale. 
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FIG. 5. Same as Fig. for N = 58 unpolarized ®He atoms. 
In the lack of TE results, we take a reference energy Eq/N 
halfway between VMCeait and LBext- Also, we do not have an 
estimate of the first excited state. Comparison of the dotted 
line with that of Fig. |^shows that the range of Et — Eo where 
the lower bound is expected to be valid shrinks as 1/A^. 



FIG. 6. Points connected by thin lines : the lower bound 
{1 — exp[— (Et — Eo)^/2(t^]}/A’ to the missing overlap per 
particle between trial function and ground state, Ct/N (with 
Eo replaced by the zero-variance extrapolation of Et)- The 
points refer to trial functions from PW to in order 

of decreasing variance. Thick lines: the approximation {1 — 
exp[—obtained using only the linear term of the fit 
of energy vs. variance. 







